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RAMAN  SIDESCATTERING  IN  LASER-PRODUCED  PLASMAS 


I.  Introduction 

Raman  Instabilities,  which  appear  in  laser-produced  plasmas  when  the 
plasma  density  is  less  than  the  quarter-critical  density,  can  lead  to  hot 
electron  production.  These  hot  electrons  could  in  turn  preheat  the  fuel  in 
laser  fusion  applications  and  are  therefore  a  significant  concern.  1-3  Of  the 
various  Raman  instabilities  in  inhomogeneous  plasmas,  sidescattering  has  the 
lowest  threshold  and  the  highest  growth  rate.  Hence,  it  is  a  natural  focus  of 
theoretical  concern. 

In  this  paper,  we  revisit  the  theory  of  this  instability.  Previous 
workers*' have  determined  the  threshold  to  be 


l“£j2  (kQL)4/3  -  1. 


where  Vq  is  the  oscillation  velocity  in  the  pump  wave  field,  kQ  is  the  wave 

number  of  the  pump  wave,  L  is  the  inhomogeneity  scale  length,  and  c  is  the 

speed  of  light.  However,  the  strictly  analytical  approach  of  Ref.  4  is  not 

sufficient  to  determine  the  accurate  numerical  value  on  the  right-hand  side  of 

Eq.  (4).  This  paper  addresses  that  issue. 

The  rest  of  the  paper  is  organized  as  follows:  Section  II  contains  the 

derivation  of  our  basic  equation.  In  Section  III  we  solve  this  equation 

numerically  to  determine  the  threshold  value.  Section  IV  contains  a  summary 

of  our  principal  results. 
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Derivation  of  the  Kenan  Sidescattering  Equation 


Our  starting  point  is  the  electron  fluid  equations  with  the  Ions  serving 
as  a  uniform  background, 
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Since  the  Instability  we  are  considering  is  parametric,  we  have  the 
approximate  relations 


“o  "  “l  +  “2 
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(3) 


where  a>  end  k  (a  -  0,1,2)  are  respectively  the  frequency  and  wave  number  of 
a  ot 

the  Incoming  pump  wave,  the  outgoing  electromagnetic  wave,  and  to  the  outgoing 
electrostatic  wave.  Equation  (3)  becomes  exact  in  the  limit  of  a  homogeneous 
plasma.  Using  the  homogeneous  limit  dispersion  relations 
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and  naglaetlng  V  /c  ,  one  finds 
e 
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(5) 


which  confirms  chat  Chis  Instability  can  take  place  only  at  densities  below 
quarter-critical . 

We  chose  our  coordinates  so  that  the  Incoming  pump  wave  Is  polarized  In 

A 

the  z-dlrectlon  and  has  Its  wave  vector  In  the  x-dlrectlon,  l.e.,  k  -  k  e  . 

X 

A  A 

We  shall  also  take  k^  -  ky  ey  -  k^  e^,  which  Is  consistent  with  sldescatterlng 

AAA 

In  an  arbitrary  direction,  so  that  k«  •  k  e  +k  e  +k  e.  We  now  write 

■*2  x  -x  y  -y  z  -z 

E  «  1  E  exp  (lk  x-w  t)  +  c.c. 

— O  X  o 


E^  -  £(x)exp  (-  Iky  y-ik^z-lw^t)  +  c.c.  (6) 

Eg  “  "  7$(x)exp  (ik^x  +  ikyy  +  lkzz  -  lwgt)  +  c.c. 

The  x-varlatlons  of  e  and  $  are  due  to  the  Inhomogeneity. 

To  proceed  we  will  use  the  following  ordering.  If  the  dispersion 


relation  for  the  a  mode  In  the  homogeneous  system  is  given  by 


D(ut  ,  k  )  -  0,  we  assume  Chat  in  Che  inhomogeneous  system  D(u>  ,  k  ),  L"'- ,  T. 

01  d  01  Ot  e 

and  nonlinearity  are  all  small,  of  order  $.  Working  to  lowest  order  in  this 
quantity  6,  we  neglect  all  products  of  these  small  quantities. 

From  Eq.  (2b)  one  finds. 


V  - 

'O 


(7) 


since  the  electromagnetic  wave  has  no  density  perturbation  and  v  is  assumed 

much  less  than  u  . 

o 

Turning  now  to  the  equation  for  the  electrostatic  wave,  we  find  from  Eq. 

(2b) 
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As  is  usually  the  case,  the  V*  TV  nonlinearity  cancels  the  VxB 

nonlinearity.  Note  also  that  we  now  use  a  subscript  2  on  the  collision 

frequency.  This  (phenomonologlcally)  allows  for  the  fact  that  the  damping  of 

the  three  different  modes  may  come  from  different  physical  processes.  For 

instance  to  account  for  Landau  damping  of  the  electron  plasma  wave,  V2  should 

be  set  equal  to  the  collision  frequency  plus  twice  the  Landau  damping  rate, 
e  eE, 

In  Eq.  (8),  Vj  ■  ,  analogous  to  Eq.  (7).  From  Eq.  (2a),  we  have 


-  lum2  -  -  Vno  V2, 


(9) 


where  n0  is  the  zero-order  background  density  which  is  x-dependent.  We  have 
also  assumed  that  nQ  and  n^  equal  zero  as  is  appropriate  for  linear 
electromagnetic  waves.  Combining  Eqs.  (8)  and  (9)  with  Poisson's  equation. 


y.  "J-  T--  3 - r^~  -3 - r-  r,  r 


Eq.  (4c),  we  conclude 
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where  v  ■  T/m.  In  obtaining  Eq .  (11),  we  replaced  to,  with  to  on  the  right- 
6  z  p 

hand  tide. 

An  analogous  derivation  yields  the  equation  for  the  sidescattered  wave 
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At  this  point,  we  must  specify  the  density  inhomogene ity.  We  shall  let 
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The  point  x  ■  0  corresponds  to  the  point  in  the  plasma  at  which  the  instability 
which  we  are  considering  occurs,  co®  is  the  plasma  frequency  at  that  point,  and  L 
is  the  density  scale  length  at  that  point.  Equations  (10)  and  (11)  then  become 

2 
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We  have  ignored  the  x-varlatlon  of  on  the  right-hand  side  of  Eq.  (13) 
consistent  with  our  ordering  of  neglecting  products  of  nonlinearities  and 
ambient  gradients. 


•V-V-V-V 
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We  further  simplify  to  the  case  of  the  scattered  wave  polarized  in  the  2- 


direction.  Since  ponderomotive  force  goes  as  ,  this  polarization  has  the 

largest  growth  rate,  at  least  in  the  homogeneous  system. 

Then  the  equations  can  be  reduced  to 
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Here  we  use  Instead  of  as  dependent  variables.  Also 
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2  17  x3x  2 


(15) 
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where  601  is  the  change  in  01  from  that  predicted  by  the  linearized  homogeneous 

2 

dispersion  relation.  We  neglect  (5oi)  ,  and  we  useoio»oi^  +oi2  to  establish 
6oij  ■  -  6oi2.  On  physical  grounds,  we  expect  that  the  threshold  will  be  lowest 
in  this  case,  and  in  the  limit  L  -►  Vj  c  >  0,  we  have  verified  that  such  is 
the  case.  Hence,  we  focus  our  attention  on  Eq.  ( 14). 

In  solving  Eq.  (14),  it  is  useful  to  take  the  Fourier  transform 

^e(x)-J  “  J  d  ^e(K)^  e3cpi  Sc*  As  We  wil1  8ee*  and  6  are  wluar« 

c 

lntegrable  and  nonsingular  along  the  real  K  axis  so  that  the  contour  may  be 
taken  as  the  real  K  axis .  Since  upon  taking  the  Fourier  transform 

Eq.  (14)  become  two  coupled  first  order  equations  in  the  K 
domain*  By  standard  means,  we  write  them  as  a  single  second  order  equation  in 


in  the  K  domain  and  eliminate  the  first  derivative  terra.  The  result  is 


+  Iy+C]  “21?+A  (19) 

d<r 

Not  only  Is  Eq.  (19)  much  simpler  Chan  Eq.  (16),  but  also  the  dimensionality 
of  the  parameter  space  has  been  substantially  reduced.  Now  y  depends  only  on 

A 

the  single  parameter  A,  whereas  in  Eq.  (16)  it  depended  on  A,  v  / c,  k  and 

e  x 

*  t 

ky.  We  will  show  that  for  parameters  of  the  NHL0  experiment,  the 
approximations  in  deriving  Eq.  (19)  are  very  well  satisifed. 


for  the  growth  rate.  It  is  interesting  to  note  that  the  limit  A  ®  does  not 
correspond  to  the  homogeneous  limit  at  x  *  0  as  one  might  naively  expect,  but 
rather  at  x  ■  (LAu^  [  iijdui  +  ioip  ^  -  i  v  ^(01°^ ]  •  It  is  possible 

to  refine^  the  estimates  in  Eqs.  (20)  and  (21)  but  an  exact  solution  of  Eqs. 


(16)  and  (19)  requires  a  numerical  approach 


We  begin  by  specifying  the  appropriate  boundary  conditions  for  E< 


Asymptotically,  as  z  goes  to  -t  «,  we  have 
<|>  -  ex  p  i  (c3/3  +  y?  )]  . 


( 


As  a  consequence  of  causality,  we  demand  that  the  solution  decay  along 
real  coasts  when  Iu(y)  >  0>  It  follows  that  we  want  the  positive  solui 
when  5  ♦  +  •  and  the  negative  solution  vtoen  5  ♦  -  •  • 

To  determine  the  eigenvalues  of  Eq.  (19),  we  calculate  the  determ 
the  solution  matrix  D(y)  using  a  finite-element  code.^  The  roots  of  tl 
determinant  D(y)  correspond  to  the  eigenvalues.  We  began  by  finding  a 
roots  of  D(y)  in  a  fairly  sizeable  box  in  the  complex  y^p  lane 
(-  6.0  <  Re  (y)  <_  6.0,  -  4.0  <_  In  (y)  <_  4.5)  with  A  set  equal  to  1.  To 
we  use  Nyquist's  theorem 


N 


<*1  P  (y  )  I 

D(y)  » 


( 


to  calculate  the  number  of  roots  N  In  the  box.  By  repeated  quadrasect 
determine  smaller  boxes  in  which  only  a  single  root  lies.  The  root  lo< 
y  is  then  estimated  using  the  relation 


v  -  1  1  LJlP.fr.>. 

Yo  2jri  *c  D(y) 


( 


The  exact  root  location  is  then  homed-in  an  using  a  secant-grid  algorli 
This  Nyquist  algorithm  is  optimized  In  a  number  of  respects  and  runs  fa 
rapidly.  (For  the  case  just  described,  the  program  tmde  roughly  35, OOt 
determinations  of  D(y)  and  took  roughly  13  minutes  c£  c.p.u.  time  to  rv 
CDC  7600.) 
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It  should  be  stressed  that  the  eigenmodes  corresponding  to  eigenvalues 
with  Im(y)  <  0  are  unphyslcal  9ince  they  correspond  to  modes  which  grow  along 
the  real  ?-axis  as  ;  ♦  ±  •,  so  that  their  Fourier  transforms  do  not  exist. 

The  eigenvalues  for  the  five  unstable  modes  which  the  routine  found  are  given 
in  Table  I  for  the  case  A  -  1.  The  number  of  significant  figures  was 
determined  by  Increasing  the  width  of  the  integration  region  and  doubling  the 
number  of  modes;  those  figures  which  did  not  change  were  deemed  significant 
and  Included  in  the  table.  The  critical  value  of  A  at  which  Im(y)  -  0,  Acrit , 
was  then  found  for  each  of  the  five  most  unstable  modes  by  following  the 
eigenvalues  parametrically  as  A  was  reduced.  The  results  are  shown  In  Table 
II,  the  overall  marginally  stability  point  is  given  by  A  ■  0.37.  We  then 

repeated  the  calculation  using  Eq.  (16)  rather  than  Eq.  (19)  with 

■*  2  2 
k  -  10,  k  «  0  corresponding  to  quarter  critical  density  and  V /c  ■  0.0013. 
x  y  6 

These  values  correspond  to  those  in  the  NRL  long  scale  length  experiment. 
Results  are  also  shown  in  Table  II.  There  is  evidently  good  agreement  between 
the  results  using  Eq.  (16)  and  the  results  using  Eq.  (19). 

In  Fig.  1  is  shown  a  plot  of  the  growth  rate  versus  A  for  the  most 
unstable  mode,  determined  from  a  numerical  solution  of  Eq.  (19).  In  physical 
units 

u°  X2 

l 2coq 6<jj  +  i Vl^  ~^T‘  "  iy(A)>  (25) 

v  1  c 

where  y(A)  is  defined  graphically  in  Fig.  1,  and  A  and  X^nt  are  defined  in 
Eqs.  (18).  A  recent  analytic  theory**  also  shows  that  the  growth  rate  depends 

on  the  single  parameter  A.  In  Ref.  8,  it  can  be  shown  that  their  expression 

1/2  -3/2 

for  growth  rate  reduces  to  y  ■  A  [1  -  0.5  A  ]  so  that  the  threshold  is 
given  by  A  *  0.63.  Thus  the  analytic  theory  predicts  a  threshold  higher  by 
almost  a  factor  of  two. 


Table  I  The  real  and  imaginary  parts  of  y  (positive  Imaginary  part  means 
growth)  for  A  ■  1  for  the  5  unstable  modes.  The  eigenvalues  are 
calculated  from  Eq.  (16)  and  also  from  Eq.  (19) 
with  kx  -  10,  ky  -  0,  3  V2/c2  -  0.00013. 


Mode  Number 

Re  y,  lay,  Eq.  (19) 

Rey,  Imy, 

Eq.  (16) 

1 

-0.4564,  0.5091 

-0.4597 

0.5168 

2 

-1.794,  0.1418 

-1.789 

0.1473 

3 

-2.835,  0.06217 

-2.827 

0.06516 

4 

-3.7114,  0.0289 

-3.7020 

0.0304 

5 

-4.49,  0010 

-4.48 

0.011 

Table  II  The  value  of  A  5  A(crit)  for  stability  for  the  5  unstable 

modes. 


Mode  Nmber 

A(crit)»  E<»*  <l9> 

A(crit)»  E<*‘  (16) 

1 

0.3  7485 

0.36762 

2 

0.6188 

0.6102 

3 

0.7567 

0.7498 

4 

0.858 

0.854 

5 

0.9 

0.9 

In  practice,  an  Important  effect  is  tranverse  convection  across  the 

finite  spot  size.  This  is  an  extremely  complicated  process.  The  theory  of 

Liu,  et  al.^  indicates  that  the  growth  rate  must  be  larger  than  the  convection 

time,  or  Im6u>  >  V^/L^  where  Lj.  is  the  transverse  scale  length.  If  VR  *  c,  (as 

it  is  at  densities  less  than  about  one-eighth  critical)  a  100  u  spot  size 

12  —1 

means  the  growth  rate  must  be  larger  than  about  3  x  10  sec  ,  a  fairly 
sizeable  requirement  on  Imdu.  However  at  the  quarter-critical  density  the 
transverse  scattered  photon  velocity  vanishes  so  that  here  the  condition 
disappears. 

From  Eq.  (25),  we  conclude  that  the  instability  threshold  is  given  by 

Im  [y(A)J  ■  w  [Vj  +  («p/Vvtl  ^Xint^c^  r4ther  than  Im  [y(A)]  -  0.  In  this 

expression,  the  v's  are  the  electron  ion  momentum  collision  frequency  plus 

twice  any  damping  rate  of  the  wave.  For  the  electron  plasma  wave  must  also 

Include  the  Landau  damping.  The  phase  velocity  of  the  wave  is  given  by 

<ii  /k».  He  have  verified  that  Landau  damping  plays  a  negligible  role  for  NHL 
po  z 

parameters  (Te  “  600  e.v.)  for  densities  larger  than  one-tenth  the  quarter- 

critical  density.  By  contrast,  the  collisions  make  a  significant  contribution 

to  the  threshold  in  density.  In  the  NRL  long  scale  length  experiment,'* 

L  •  140y,  Z  **  5 ,  T  *00  eV.  Thus  at  the  quarter-critical  density,  taking 

into  account  collision  damping,  the  threshold  for  Raman  sidescatter  is  A  - 

0.46  or  I  ■  1.6  x  10**w/cn2  using  the  fact  that  V  •  25.6  I^2(w/cm2)  X  (y). 

os  o 

IV.  Summary 

In  this  paper  we  have  revisited  the  theory  of  Raman  sidescattering  in 
laser  produced  plasmas.  We  have  shown  that  to  a  very  good  approximation  the 
growth  rate  and  thresholds  can  be  reduced  to  examining  a  single  parameter 


system.  We  find  that  at  quarter-critical  density,  the  threshold  for 


stimulated  Raman  sldescatter  In  the  NHL  long  scale  length  experiment  Is  given 
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roughly  by  1  >  1.6  x  10  m/cm  • 
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